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В данной работе описывается простой и эффективный алгоритм, реали- 
зующий замкнутый набор булевых операций {объединение, пересечение, 
разность, симметрическая разность} над множествами гранично задан- 
ных многоугольников (Боипаагу гергезете4 ро[уеоп5) на плоскости. Замк- 
нутость набора операций подразумевает корректную обработку много- 
угольников с самокасаниями (или кратными вершинами), т.е. результи- 
рующие многоугольники, полученные с помощью предлагаемого алгорит- 
ма, удовлетворяют предъявляемым требованиям к входным данным алго- 
ритма. Асимптотическая оценка времени выполнения алгоритма не превы- 
шает О(и [о2*п + +2021), где п — общее число вершин многоугольни- 
ков-операндов, К — число точек пересечения их рёбер, 2 — общее число 
контуров многоугольников-операндов. Таким образом, предложенный ал- 
горитм по эффективности не уступает большинству алгоритмов, известных 
на сегодняшний день. 
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ВВЕДЕНИЕ 


Каждый простой двумерный многоугольник ограничивает некоторую 
область на плоскости. Под выполнением булевой операции 
ор Е {объединение, пересечение, разность, симметрическая разность} над 
такими многоугольниками понимается выполнение операции ор над опи- 
сываемыми ими областями и получение множества многоугольников, опи- 
сывающего полученную область. Основная проблема заключается в том, 
что результатом некоторой операции даже над двумя простыми выпуклыми 
многоугольниками в общем случае может являться множество невыпуклых 
многоугольников с возможными отверстиями и самокасаниями (например, 
А - В на рис. 1). 

К сожалению, традиционные алгоритмы 
ограничиваются тем, что разбивают 
сложные многоугольники результата на 
простые, тем самым исключая возмож- 
ность повторного выполнения операций 
над результатом. Для задач отсечения 
по экрану и выяснения критических 
участков интегральных схем этого дос- 
таточно, но в системах  автоматизи- 
рованного проектирования и управле- 
ния географическими базами данных 
требуется реализация замкнутого набора операций над множествами мно- 
гоугольников. Это означает, что результат, полученный с помощью алго- 
ритма, должен удовлетворять его же предусловию. Также желательно, что- 
бы алгоритм допускал отверстия и неограниченную кратность вершин. 
Проведём краткий обзор существующих алгоритмов. 

Вейлер (У/еПег) и Азертон (АФеноп) [5, 12, 13] предложили следующий 
алгоритм. Ребра исходных многоугольников дублируются, и между ними 
вычисляются пересечения. Затем производится перераспределение ссылок 
в вершинах пересекающихся ребер с целью выделения минимальных кон- 
туров. Для этого используется довольно сложный набор правил для различ- 
ных случаев пересечения рёбер. На основе полученных контуров создаются 
результирующие многоугольники для операций объединения, пересечения 
и вычитания. Перечислим основные недостатки этого подхода: 

— результат булевой операции не всегда состоит из минимальных кон- 
туров (см. 0, результат операции вычисления симметрической разности 


АВ); 


Рис. 1 


— не освещена проблема описания многоугольников с самокасаниями, 
которые могут быть получены в результате работы алгоритма; 

— результат операции может не удовлетворять требованиям на входные 
данные, причем расширение алгоритма для обработки кратных пересечений 
крайне затруднительно. 

Шутте (Зевиаве) [9, 19] модифицировал алгоритм Вейлера более четким 
разделением этапов и новым алгоритмом маркировки ребер (еде 1абепз), 
однако его алгоритм обладает более существенными ограничениями на 
входные данные. 

Итак, приведенные выше алгоритмы реализуют набор операций, не яв- 
ляющийся замкнутым, так как в результате их выполнения могут получать- 
ся контуры с самокасаниями, что неприемлемо в качестве исходных дан- 
ных для обоих алгоритмов, и с отверстиями, которые не допускаются алго- 
ритмом Шутте. 

Сравнительно недавно Гютинг (Счейп?) и Шнайдер (Зсбпе!4ег) [6 - 8] 
разработали КОЗЕ-алгебру (ВОБи${ Брайа1 Ежепзюп а{оебга), реализую- 
щую замкнутый набор операций над двумерными объектами, в том числе 
регионами (гег1оп$). Предложенное ими описание также имеет ряд недос- 
татков: 

— координаты вершин КОЗЕ-объектов принадлежат дискретному мно- 
жеству, в отличие от общепринятой в вычислительной геометрии модели 
вещественнозначной РАМ [13]; 

— сложный и громоздкий в реализации алгоритм связывания сегмен- 
тов (Ва-зеотеп($) региона в корректные внешние и внутренние циклы (В- 
сус[ез), приводящий к снижению эффективности процедур обработки ре- 
ГиОНОвВ. 

В данной работе описывается простой и эффективный алгоритм, кото- 
рый свободен от перечисленных выше недостатков и реализует замкнутый 
набор операций {объединение, пересечение, разность, симметрическая 
разность} над множествами многоугольников на плоскости. В разд. 0 да- 
ются основные определения, в разд. 0 описывается собственно алгоритм, в 
последующих разделах анализируется его эффективность и обсуждаются 
некоторые детали реализации. 


1. ОПРЕДЕЛЕНИЯ И ОБОЗНАЧЕНИЯ 


Далее ребром называется направленный отрезок с началом в точке а и 
концом в точке Б, не совпадающей с а, котрый обозначается как Ё(а, Б). 
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При этом будем говорить, что для точки В ребро Ё является входящим, а 
для точки а — выходящим. Множество точек х, удовлетворяющих услови- 


ям 0 < /(Бах) < ф и 2п-ф < И(абх) <П2лт, назовем левой Ф-окрестностью 
ребра Е'. Множество точек х, удовлетворяющих условиям 0 <1(аБх) < 


и 2^-Ф < И(Бах) <О2лт, назовем правой Ф-окрестностью ребра Е (см. рис. 
2а). 

Определение 1. Контуром называется упорядоченное множество п ре- 
бер С = {ЁЕь, Еф, ..., Е п, п>3 такое, что УТЕ {0...п-1} Её: =Е(уь м) и 
Е=Е(Уь а 

Точку у, для которой ребра Е; | и Е; являются соответственно входящим 
и выходящим, будем называть {-й вершиной контура С, ребра ЕЁ; 1 и Ё; — 
соседними или смежными. Порядок обхода ребер контура С = {Еь, Ё1, ..., 
Е 1} с увеличением индекса (1-—> 1+ 1) будем называть прямым направле- 
нием обхода, а противоположный ему — обратным. 

Заметим, что одной точке плоскости может соответствовать несколько 
вершин контура, такие вершины будем называть кратными. Углом входя- 
щего в вершину у ребра Е(а,») будем называть направленный угол между 
векторами е! (орт абсцисс) и Е(у,а). Углом выходящего из вершины у ребра 
Е(у,Б), будем называть направленный угол между векторами е! и Е(»,Ь). 

Рассмотрим пару соседних ребер Е; | и Ё; контура С: 0 <1< и — 1. Мно- 
жество точек х, удовлетворяющих условиям 0 <П1/(ун их) < Хил) и 


О<р-уъ| <&, назовем внутренней 8-окрестностью вершины у; и будем 
обозначать 5'(”). Множество точек х, удовлетворяющих условиям 0 < 
Ионы х) < Ими у) ид < К- у! < 8 назовем внешней 8-окрестностью 
вершины ъ, и будем обозначать & (у) (см. рис. 26). 

Рассмотрим евклидово пространство размерности 2, наделенное декар- 
товыми координатами с ориентирующим репером ег, е»›, т. е. плоскость с 
традиционным “правым” репером. 


Определение 2. Областью называется ограниченное полигональное от- 
крытое связное множество на евклидовой плоскости, определенное с точ- 
ностью до множества меры нуль. 


| Г (АВС) обозначает направленный угол между векторами ВА и ВС, лежащий в интервале 
[0,2п). 

? Здесь и далее при вычислении индексов вершин и ребер используется сложение и вычи- 
тание по модулю п, где и — число ребер рассматриваемого в данный момент контура. 
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Левая (р-окрестность ребра Е(а, 5) 


ф- 


Правая ф-окрестность ребра Е(а, 5) 


6) У 


У ы = 
ВА. 


Рис. 2 

Под полигональностью понимается то, что граница области состоит из 
замкнутых ломаных. Определение области с точностью до множества меры 
нуль позволяет рассматривать области без “разрезов” и “дыр” в точке. За- 
метим, что область может быть как одно-, так и многосвязной. 

Рассмотрим границу АА некоторой области А. Если ориентация АА со- 
гласована с ориентацией А, мы приходим к соответствию между границей 
области и набором контуров на плоскости. Для того, чтобы оно было вза- 
имно-однозначным, добавим к традиционному определению два ограниче- 
ния. Итак, рассмотрим контур С, содержащий и ребер. 

Определение 3. Контур С называется ограничивающим контуром об- 
ласти А, если выполняются следующие условия: 

П СеаА; 

2) обход С в прямом направлении является положительным обходом, то 
есть обходом, при котором область А остается слева; 

3) УЕе {0..п-13 &>0:У хеё) хед; 
УТЕ {0...п- У &>0:3 хеЕБ (у) х®А. 


Определение 4 //олигоном называется совокупность всех ограничи- 
вающих контуров некоторой области. 

Внешним контуром области А назовем ограничивающий ее контур, по- 
ложительный обход которого производится против часовой стрелки. Внут- 
ренним контуром области А назовем ограничивающий ее контур, положи- 
тельный обход которого производится по часовой стрелке. Очевидно, что в 
произвольном полигоне может быть один и только один внешний контур. 


Определение 5. Регионом называется набор полигонов, описывающий 
множество непересекающихся областей. 


Рис. 3 
Таким образом, установлено взаимно-однозначное соответствие между 
множеством непересекающихся областей и множеством ограничивающих 
их контуров. На рис. 3 можно увидеть удовлетворяющие данным выше оп- 
ределениям описания областей (закрашенные) регионами (пунктирными 
линиями показаны направление обхода и связность ограничивающих кон- 


туров). 


2. ОПИСАНИЕ АЛГОРИТМА 


Входные данные алгоритма: 2 региона А и В, операция ор = 
{ объединение, пересечение, разность, симметрическая разность}. 
Выходные данные алгоритма: регион К = А ор В. 
Алгоритм состоит из четырех этапов. 
1. Обработка пересечений ребер. 
2. Маркировка контуров и ребер. 
3. Получение результирующих контуров. 
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4. Создание результирующего региона. 

Рассмотрим в качестве примера регионы А и В, изображенные на рис. 4. 
С учетом наших определений регион А состоит из внешнего и внутреннего 
контуров, а регион В — из внешнего контура, касающегося самого себя в 
точке. 


Рис. 4 


2.1. Обработка пересечений ребер 


Найдем все пары не смежных пересекающихся ребер регионов А и В. 
Если пересечение ребер не точка, а отрезок, будем считать концы этого 
отрезка двумя точками пересечения. Если точки пересечения не совпадают 
с концевыми точками ребер, в исходные регионы добавляются новые вер- 
шины с координатами соответствующих точек. Отметим, что добавление 
точек пересечения не изменяет области, описываемые регионами А и В. 
Далее везде под регионами А и В имеются в виду исходные регионы А и В 
со всеми добавленными точками пересечения. Вершины, соответствующие 
точкам пересечения ребер, назовем вершинами-пересечениями (отмечены 
жирными точками на рис. 5). Отметим, что одной геометрической точке 
пересечения ребер соответствуют как минимум две вершины-пересечения. 

Рассмотрим вершину-пересечение у одного из регионов. Если за Е"(ъ) 
обозначим входящее в у ребро, а за Ё (у) — выходящее, тогда обработка 
вершины-пересечения выполняется следующим образом. 

1. Создаются два дескриптора вершин-пересечений О*(у) и )(», отве- 
чающие соответственно Е*(у) и Е (ъ). В дескрипторы помещается информа- 
ция об углах их ребер, а также о том, являются ли эти ребра входящими или 
выходящими относительно у. 
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2. Между у и дескрипторами О*(у) и О’ (») устанавливаются двусторон- 
ние связи, которые позволяют определить соответствующую некоторой 
вершине-пересечению пару дескрипторов и наоборот. 

3. 2*(»›) и >) помещаются в список связности [(х) точки пересечения 
х. Данный список является кольцевым списком дескрипторов вершин- 
пересечений, соответствующих точке х, отсортированным по углам ребер. 
В дальнейшем список Г(х) будет использоваться для поиска ребер, бли- 
жайших к данному в направлениях по часовой и против часовой стрелки. 
Порядок в списке связности дескрипторов с одинаковыми углами ребер не 
специфицирован. 


Рис. 5 
На рис. 6 приведен пример списка связности точки х, выделенной на 
рис. 5. Ей соответствуют три вершины-пересечения — а, В» и Во. 


т: В В В.В В В. 


Рис. 6 
Далее будем говорить, что вершина-пересечение у присутствует в неко- 
тором списке связности, если в нем находятся ее дескрипторы О*(») и Б`(ъ). 


И 


Так как у может присутствовать только в списке связности геометрически 
соответствующей ей точки, можно употреблять понятие списка связности 
вершины-пересечения. 


2.2. Маркировка контуров и ребер 


Рассмотрим некоторый ограничивающий контур С одного из регионов 
Аи В. Обозначим за М регион, к которому не принадлежит С. 

На предыдущем этапе все точки пересечения были добавлены в исход- 
ные регионы. Отсюда, с учетом определений 3, 4 и 5, с очевидностью выте- 
кает следующее утверждение. 


Утверждение 1.После выполнения первого этапа алгоритма: 

1) ребра регионов А и В не имеют никаких общих точек друг с другом, 
кроме концевых; 

2) произвольное ребро контура С может либо геометрически совпа- 
дать с одним из ребер М (такие ребра будем называть сопряженными), либо 
находиться целиком (за исключением, может быть, концевых точек) внутри 
или вне области, описываемой М; 

3) если контур С не содержит вершин-пересечений, то он лежит цели- 
ком внутри или вне области, описываемой М. 

Меткой ребра Е, принадлежащего контуру С, будем называть его атри- 
бут, имеющий значение из множества {5НАКЕР!, 5НАКЕР2, ПМ5ШЕ, 
ОИТУ5ШЕ\ и отражающий его геометрическое положение относительно 
региона М следующим образом: ИМ5/Ш)Е — ребро Е (может быть, за исклю- 
чением его концевых точек) лежит внутри М, ООТУШЕ — ребро Е лежит 
вне М, ЗНАКЕР! — ребро Е и сопряженное ему ребро из М сонаправлены, 
5НАКЕР?2 — ребро Е и сопряженное ему ребро из М направлены в проти- 
воположные стороны. 

Меткой контура С мы будем называть его атрибут, имеющий значение 
из множества {/5ЕСТЕД, Г[М5ШЕ, ОПТ5ШЕ} и отражающий его геометри- 
ческое положение относительно М следующим образом: /5ЕСТЕР — со- 
держит вершины-пересечения, М5 1)Е — С лежит внутри М, ОЧТУШЕ — 
С лежит вне М. 

Корректность определений меток контура и ребра следует из утвержде- 
ния [. 

Приведем алгоритм маркировки контура С и его ребер. 

1. Сне содержит вершин-пересечений. 
Если С лежит внутри М, то С помечается как /М5УШЛЕ, в противном слу- 
чае — как ООТ5Е. Маркировки ребер С не производится. 
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2. С содержит вершины-пересечения. 

С помечается как [5ЕСТЕР и производится последовательная марки- 
ровка всех его ребер Е, ге {0...м-1}, где и — число ребер в С. Пусть Еха, 
Ь) — маркируемое ребро. 

2.1 Е; не содержит вершин-пересечений. 
1) 1+0. Значение метки копируется у ребра Е; 1. 
2) 1=0. Если а лежит внутри М, то Ё; помечается как //М5/ДЕ, в про- 
тивном случае Е; помечается как ОПТ5У/Е. 
2.2 Е; содержит вершины-пересечения. 
1) Зребро Ес, 4: Е е Мла=сльЬ=а. Е; помечается как 5НАКЕР/. 
2) Зребро Е(с, а: Е е Мла=алЬ=с. Е, помечается как 5НАКЕО2. 
3) В списке связности вершин-пересечений не присутствует вершин 
М. 
Такое возможно, если контур С не пересекается с М и касается са- 
мого себя. Метка ставится так же, как и в случае 2.1. 

4) В списке связности вершин-пересечений присутствуют вершины 
М. 
Для каждой из таких вершин у^, определяется, лежит ли Ё; внутри 
И имил). Если Е, попало во внутренность хотя бы одного из 
таких “маркирующих” углов, то оно лежит внутри М и помечается 
как ИМЕ. В противном случае Е; помечается как О(Т5ШЕ(см. 
рис. 7). 


Ребро помечается как 
ОПТ5ШЕ, т.к. оно лежит 
вне “маркирующих” углов 


Ребро помечается как 
М$5ШТЕ, т.к. оно лежит 
внутри “маркирующего” 
угла 


м ----- — 
Маркируемые ребра Ребра региона М 
Рис. 7 
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2.3. Получение результирующих контуров 


2.3.1. Предварительные рассуждения 


Одна из ключевых идей данного алгоритма заключается в том, что ог- 
раничивающие контуры региона РА (далее результирующие контуры) соби- 
раются на основе меток контуров и ребер, а не координат их вершин. Так 
как одно и то же ребро не может входить в регион более одного раза, нам 
нужно найти правила, по которым ребра и контуры включаются в резуль- 
тат. 

На рис. 8 приведен результат маркировки регионов А и В нашего при- 
мера. 


А 
ОЦТЗШЕ 
—4А-^ 


М$5ШЕ 
— 


5НАКЕР2 


— 


5НАКЕБ1 


—ъ 


Рис. 8 


2.3.1.1. Условия включения ребер. 
Рассмотрим некоторое ребро Е, принадлежащее А или В. Выберем наи- 


большее Ф > 0 такое, что ни левая, ни правая Ф-окрестности ребра Е не 

будут содержать никаких вершин и ребер регионов А и В. Назовем эти ок- 

рестности левой и правой минимальными окрестностями ребра Е и обо- 
+ 2 

значим их соответственно Ф (Е) и Ф (ЕЁ). По выбору Ф все точки мини- 
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мальной окрестности ребра одновременно принадлежат или не принадле- 
жат области, описываемой регионом А или В. 

Принадлежность некоторого ребра результату можно установить по 
принадлежности результирующей области минимальных окрестностей реб- 
ра следующим образом: 


+ - 2 
1. Если Ф (ЕЁ) входит в результат, а Ф (Е) — нет, Е войдет в результат 
с исходной ориентацией. 


2. Если Ф (Е) входит в результат, а ФЕ) — нет, Е войдет в результат 
с обратной ориентацией. 

3. Если обе минимальные окрестности Ё не входят в результат, Е не 
войдет в результат. 

4. Если обе минимальные окрестности Ё входят в результат, Ё не вой- 
дет в результат (иначе возникает противоречие пункту 4 определения 3). 

Заметим, что из пары сопряженных ребер в результат можно включить 
только одно с нужной ориентацией. 

В табл. 1 рассмотрены условия включения для всех возможных случаев 
меток ребер обоих регионов. 


Таблица 1 


Случай 1. 
Е принадлежит региону А и имеет мет- 


ку ОПТУШДЕ; ФЕ) принадлежит региону 


А, Ф (Е) не принадлежит ни одному из 
регионов; Ё войдет в результат операций 
объединение, разность и симметрическая 
разность с исходной ориентацией. 
Случай 2. 
Е принадлежит региону В и имеет мет- 


ку ООТУШЕ; ФЕ) принадлежит региону 


В, Ф (Е) не принадлежит ни одному из 
регионов; Ё войдет в результат операций 
объединение и симметрическая разность 
с исходной ориентацией. 
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Случай 3. 
Е принадлежит региону А и имеет мет- 


ку ПМ5ШЕ; ФЕ) принадлежит обоим ре- 


гионам, ф< (Е) принадлежит региону В; Е 
войдет в результат операции пересечение 
с исходной ориентацией, а в результат 
операции симметрическая разность — с 
обратной. 

Случай 4. 

Е принадлежит региону В и имеет мет- 


ку ПМ5ШЕ; ФЕ) принадлежит обоим ре- 


гионам, ф< (Е) принадлежит региону 4; Ё 
войдет в результат операции пересечение 
с исходной ориентацией, а в результат 
операций разность и симметрическая 
разность — с обратной. 

Случай 5. 

Пара двух геометрически совпадаю- 
щих и направленных в разные стороны 
ребер ЕЕ АиЕЕ В с метками 5НАКЕР2; 


Фф’(Е) принадлежит только региону А, а 


Фф (Е) принадлежит только региону В. 
Одно из ребер войдет в результат опера- 
ции разность с направлением, совпадаю- 
щим с направлением ребра Е. 

Случай 6. 

Пара двух геометрически совпадаю- 
щих и направленных в одну сторону ре- 
бер регионов А и В с метками 5НАКЕБ/. 
Левая минимальная окрестность ребер 
принадлежит обоим регионам, а правая — 
ни одному из них. Одно из ребер войдет в 
результат операций объединение и пересе- 
чение с исходным направлением. 


2.3.1.2. Процедура перехода в вершине-пересечении 
Если в процессе сборки ребер результата мы пришли в вершину- 
пересечение у, возникает проблема выбора дальнейшего направления обхо- 
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да. Пусть х — точка плоскости, соответствующая у. Рассмотрим все верши- 
ны-пересечения, геометрически совпадающие с х. Пусть их число равно т, 
т > 2. Тогда число их дескрипторов равно 2т (см. 0). С помощью Г(х) по- 
строим множество ребер РЁ’, /е{0...2т-—1} такое, что Ру = Е*(у) и дескрип- 
тор, соответствующий ребру Ё.-! является ближайшим в направлении по 
часовой стрелки к дескриптору, соответствующему ребру Р;. 


Выберем наибольшее $ > 0 такое, что ни одна из внутренних или внеш- 


них 5-окрестностей рассматриваемых вершин не будут содержать никаких 
вершин и ребер, отличных от Е. Далее рассмотрим сектора, на которые 


геометрическая 5-окрестность точки х разбивается ребрами РЁ.. Назовем 
такие сектора минимальными. Очевидно, что все точки минимального сек- 
тора (за исключением, может быть, его границ) одновременно принадлежат 
или не принадлежат области, описываемой регионом А или В. Если точки 
соседних минимальных секторов должны войти в результат, то по опреде- 
лению 3 их должен ограничивать один и тот же контур результирующего 
региона. 

Таким образом, нужно, начиная с Ру и перебирая 2; в направлении воз- 
растания /, найти ребро, удовлетворяющее условию включения для выпол- 
няемой операции. Оно и будет искомым ребром для продолжения обхода. 
Отметим, что для поиска не требуется информации о координатах вершин, 
так как дескрипторы отсортированы в списке связности по углам соответ- 
ствующих ребер. 


2.3.2. Алгоритм сборки результирующих контуров 


Так как предыдущие рассуждения справедливы для всех операций, 
сборка результата осуществляется для всех операций единообразно. На 
данном этапе последовательно рассматривается каждый контур из регионов 
Аи В. Обозначим рассматриваемый контур за С. 

Метку контура или ребра Х будем обозначать Х.ЁЕ1аэ5. Далее обход ре- 
бер контура в прямом и обратном направлениях (см. определение 1) обо- 
значим как КОК ТАКО и ВАСКУ’АЮР соответственно. Приведем алгоритм 
сборки результирующих контуров. 

1. Если С. Нае5 = 5ЕСТЕР, то контур включается в результат на осно- 
ве своей метки с нужным направлением. В табл. 2 приведены условия вклю- 
чения для контура С в зависимости от выполняемой операции ор. 


И 


Таблица 2 
ор Условие включения С в результат Направление 


Ю (С.Е1ав5 = МЕ) РОВТАЕО 
№ (С.Е/ав5 = ООТУЮЕ) РОКТАЕЬ 


№ (С.Наз5 = ОИТУ$ШЕ) ^ (СЕА) ЕОЕТАКО 
(С.Наз5 = МУШЕ) ^ (СЕ В) ВАСКТАКЬ 
= (С.Наз5 = ОИТУШЕ) ЕОЕТАКО 
(С.Наэ5 = МУШДЕ) ВАСКИАКО 
2. Если С.Е/азз = [5ЕСТЕО, нужно совершить обход всех ребер С, что- 
бы найти ребра, с которых начнется сборка результирующих контуров. Ес- 
ли п — число вершин контура С, то алгоритм поиска подходящих ребер 
выглядит следующим образом: 
Гогт := 0 10 п-1 40 
Бест 
И (ЕазеВще(Е, 4’) апа (Е; не помечено как пройденное)) 
И (4: = ЕОВУАВО) 
г := СоПесК(», а); 


е15е 
г := СоПесу.1, ай); 
Включить г в набор результирующих контуров; 
епа; 
где РапеНоп ЕдгеВше(Е:ребро; уаг даи:(ЕОВУАВО, ВАСКУАВО)): 
Бооеап — функция-параметр, соответствующая условию включения ребра 
в результат выполняемой операции (см. табл. 3). 
Таблица 3 
ор Условие включения Е в результат Направление 


г | (Е.Е/ав5 = МУШЕ) МС (Е.Е а = ЗНАКЕБ) КОКТАКО 
кл | (Е.Р 1985 = ООТУШЕ) МС(Е.Е1ая5 = 5НАКЕОГ) КОКТАКО 


— | (Е Наву = ООТЯШЮЕ) РОКТАЕО 
у (Е.Е1авз = 5НАКЕР?)) л (ЕЕ А) 
((Е.Назз = МУШЕ) ВАСКТАКО 
и (Е.Е/ае5 = 5НАКЕР?)) л (Е Е В) 

Ф | (ЕРа = ООТЯЮЕ) РОВТАКР 


(Е.Н1аэ5 = М№5ШЕ) ВАСКУАКО 
Процедура сборки результирующего контура 7, начиная с вершины у, 
выглядит следующим образом: 
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Гапесйоп СоПесКу:вершина; 9г:(ГЕОВУАВО, ВАСКУ/АВО)):контур; 
Безт 
Создать пустой контур г; 
гереае 
Включить увг; 
И (4 = ЕОВУУАВО) 
Е := выходящее из У ребро; 
е15е 
Е := входящее в У ребро; 
Пометить Е как пройденное; 
И ((Е.Р/аэ5 = 5НАКЕОТ) ог (Е.Е 1аэ5 = ЗНАКЕО2)) 
пометить сопряженное Е ребро как пройденное; 
у := вершина, следующая относительно У в направлении а; 
И (у является вершиной-пересечением) еп 
Литр(у, 41); 
ипё (текущее ребро помечено как пройденное); 
гебиги г; 
епа; 


ргоседиге Латр(уаг у:вершина; уаг 4: (ЕОВУАВО, ВАСКУУАВО)); 
Бест 
Н (4: = ГОВУУАВО) еп 
Я := ргеу(0“(у)); 
ее 
4 := ргеу(О `(у)); 
{За ргеу(О) мы обозначили дескриптор, ближайший по часовой стрел- 
кеко} 
Гоипа := ЕКАГЗЕ; 
гереаё 
е := ребро, соответствующее @; 
Н ((е не помечено как пройденное) ап ЕдзеВше(е, пеуудй`)) еп 
Берт 
И ((е входящее) ап4 (пеуудй: = ВАСКУУАВО) ог 
(е выходящее) ап4 (пеуудй: = ЕОВУУАВО)) еп 
Бест 
г := пеуд; 
у := вершина-пересечение, соответствующая 4; 
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Гоипа := ТВОЕ; 


 := ргеу(а); 
ип (0ип9); 
епа; 


На рис. 9 приведены результирующие контуры различных операций над 
регионами А и В примера. 

Результат операции объединение состоит из одного внешнего и двух 
внутренних контуров, разность и пересечение — из двух внешних конту- 
ров, симметрическая разность — из трех внешних контуров. 


2.4. Создание результирующего региона 


Итак, мы получили набор ограничивающих контуров А. В отличие от 
многих других алгоритмов, задача правильного формирования А очень про- 
ста, так как полученные на предыдущем этапе внешние и внутренние кон- 
туры имеют правильную ориентацию. 

Для каждого внешнего контура создается отдельный полигон, внешним 
контуром которого он и будет являться. Внутренние контуры на предыду- 
щем этапе сохраняются во временном списке, а затем помещаются в поли- 
гоны результата. Найти нужный полигон несложно, для этого достаточно 
определить минимальный из внешних контуров, содержащих данный внут- 
ренний контур. Регион А будет состоять из множества полученных таким 
образом полигонов. 
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3. АНАЛИЗ ЭФФЕКТИВНОСТИ АЛГОРИТМА 


Проанализируем эффективность описанного выше алгоритма. Пусть ис- 
ходные регионы А и В содержат 
всего п вершин и 2 контуров. 

На этапе 0 максимальное число 
новых вершин А составляет О(и?) 
(см. рис. 10). Время выполнения 
данного этапа определяется исполь- 
зуемым алгоритмом нахождения 
пересечений. Вычислительная сло- 
жность наиболее простого алгорит- 
ма попарного перебора ребер со- 
ставляет О(и? + К). Существует ряд Рис. 10 
более эффективных алгоритмов 
[1 -3], время выполнения которых составляет О(и [ое п- ^). В работе Фин- 
ке (ЕшКе) и Хинрикса (Ншисвз) [4] описывается алгоритм, который можно 
использовать для нахождения пересечения двух многосвязных подразбие- 
ний плоскости (раф  соппеме4 р|апаг  за691\15100$) за время 
О(п 10о5*п+ К+210е п) (пусть п>0, [0°®и=и, [ое®и = [о2(ю®1и), тогда 
102*п — это наибольшее целое число 5, такое что /ое®и > 1). 

Этап 0 выполняется за время О(и +=р(и)), где р(п) — время, необходи- 
мое для выяснения, находится ли некоторая точка внутри региона, состоя- 
щего из и вершин. Без использования дополнительных структур данных 
р) = О(т). Для ускорения выполнения таких запросов можно использо- 
вать структуру, полученную при выполнении  трапезоидации 
(Барего14айоп) регионов, которую можно построить за О(и [оэ*и + 2 [08 п) 
11]. С помощью этого приема можно достичь оптимального времени 
Р(п) = О(о® п). 

Так как на этапе 0 каждое ребро проходится не более двух раз, он вы- 
полняется за время О(и + К). Время выполнения этапа 0 не превышает 
О(= рп + К). 

За счет использования минимальных описывающих прямоугольников, 
которые вычисляются за время О(и), можно достичь дополнительного 
уменьшения вычислительных затрат на всех этапах алгоритма. 

Таким образом, время выполнения алгоритма не превышает 
О(п [02*п+ К + 202 п), и затраты на память, очевидно, составляют О(и + ^). 
Эти результаты на сегодняшний день являются оптимальными как с теоре- 
тической, так и с практической точки зрения. 
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ЗАКЛЮЧЕНИЕ 


Перечислим основные результаты данной работы. 

® Даны четкие определения, позволяющие однозначно описать полиго- 
нальные области на плоскости с отверстиями и самокасаниями. 

®Разработан метод обработки кратных вершин-пересечений, позволив- 
ший решить проблему вырожденности (4есепегасу) точек, характерную для 
традиционных алгоритмов. 

®Введен набор ясных и интуитивно понятных (обращаем внимание на 
симметрию условий включения в результат контуров и ребер) правил для 
сборки результирующих контуров, допускающих эффективную и компакт- 
ную реализацию алгоритма. 

Главным достижением данной работы является замкнутость получен- 
ного набора операций. 

Авторская реализация алгоритма занимает около 850 строк на языке С, 
причем в случае больших полигональных моделей реальное время выпол- 
нения несколько меньше теоретически предсказанного ввиду использова- 
ния различных оптимизирующих техник (см. разд. 0). 

В заключение хочется поблагодарить Ральфа Гютинга, Раймунда Зейде- 
ля и Джека Сноуинка за предоставленные материалы, и Кламера Шутте, с 
готовностью обсуждавшего с авторами плюсы и минусы своего алгоритма. 
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